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摘要 : 随 着 射电 天 文 的 发 展 ， 大 口径 天 线 的 角 分 状 率 越 来 越 高 ， 对 天 线 指向 精度 的 要 求 
高 达 1" 量 级 。 由 于 天 线 伺服 驱动 力矩 的 脉动 造成 低频 谐振 ， 传统 PI(D) 调 节 器 难以 满足 更 
高 指向 跟踪 的 需求 。 结 合 射电 望远镜 的 工 况 、 伺 服 控制 目的 及 控制 器 设计 原则 ， 采 用 一 种 基 
于 线性 自 抗 扰 控制 算法 ， 并 应 用 于 某 天 线 位 置 环 仿真 。 分 析 结 果 表 明 ， 该 算法 不 仅 有 良好 的 
动态 性 能 和 控制 精度 ， 而 且 有 较 强 的 重 棒 性 。 
关键 词 ， 指 向 精度 ; 低频 谐振 ; 伺服 控制 ; 线性 自 抗 扰 控 制 ; 鲁 棒 性 
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全 可 动 射 电 望 远 镜 主要 应 用 于 射电 天 文 观测 和 深 空 探测 ， 准 确 指 向 并 稳定 跟踪 目标 天 体 是 射电 望 
远 镜 的 基本 功能 。 随 着 射电 天 文 研究 的 深入 和 天 线 制 造 技 术 的 发 展 ， 全 可 动 射电 望远镜 的 口径 越 来 越 
bm 大 ， 天 线 重 量 随 之 增加 ， 但 天 线 结构 刚度 下 降 ， 进 而 使 得 天 线 固 有 频率 降低 ， 对 天 线 控制 精度 带 来 巨 
rr 大 的 挑战 。 

根据 天 线 接收 信号 功率 不 低 于 50% 及 分 辩 要 求 ， 定 义 天 线 的 角 分 辩 率 为 天 线 半 功 率 束 宽 (HPBW ) 。 
~N 为 了 使 天 线 指向 引起 的 增益 误差 小 于 3%， 要 求 天 线 电 轴 指 向 精度 小 于 半 功 率 束 宽 的 1/10。 在 跟踪 模 
= 式 下 ， 为 使 望远镜 电 轴 保持 在 指向 精度 内 ， 要 求 望远镜 伺服 系统 具有 优越 的 控制 精度 。 

目前 ， 观 测 100 GHz 以 上 波段 时 ， 美 国 100 m 绿 岸 望远镜 ( Green Bank Telescope, GBT) 要 求 平均 
风速 小 于 1 m/s 时 的 指向 精度 优 于 1"， 平均 风速 小 于 4m/s 时 的 指向 精度 优 于 1.5”; 国内 65 m REE 
远 镜 要 求 风速 小 于 4 m/s 时 的 指向 精度 优 于 3"。 拟 在 新 疆 奇 台 建 设 的 110 m 全 可 动 射电 望远镜 (QiTai 
radio Telescope，QTT)' 要 求 指向 精度 优 于 1. 5", 满足 7 mm 甚至 3 mm 波段 观测 对 天 线 指向 精度 的 要 
求 ， 无 疑 对 伺服 控制 设计 提出 极 高 的 挑战 。 


O ”1 射电 望远镜 天 线 指向 的 影响 因素 


射电 望远镜 天 线 电 轴 通过 机 械 轴 带动 ， 对 于 轮 轨 式 天 线 ， 机 械 轴 靠 天 线 方位 电机 和 俯仰 电机 的 运 
动 实现 。 理 想 情况 下 ， 天 线 电 轴 与 机 械 轴 重 合 。 由 于 加 工 精度 、 安 装 工艺 ”、 天 线 自重 、 工 作 环境 
温度 变化 、 雨 雪 等 因素 使 得 天 线 结构 形变 ， 造 成 电 轴 与 机 械 轴 偏差 ， 且 不 同 俯仰 角 和 大 气 折射 率 '” 
都 影响 天 线 的 电 轴 指向 ， 这 些 确定 性 载荷 可 通过 误差 修正 函数 修正 。 影 响 电 轴 指向 的 随机 而 复杂 的 因 
素 主 要 是 天 线 运 行 过程 中 受到 的 齿 隙 非 线性 、 天 线 柔 性 、 系 统 噪 声 、 系 统 参 数 的 摄 动 和 阵风 载荷 。 研 
究 发 现 ， 天 线 在 低速 运行 时 ， 由 于 天 线 机 械 部 件 刚度 不 一 致 以 及 路 合 齿轮 间 存 在 的 库伦 摩 控 ， 电 机 驱 
动 天 线 架 的 力矩 发 生 脉 动 变化 ， 而 阵风 作用 的 实质 是 给 天 线 系统 额外 施加 了 一 个 低频 脉动 力矩 ， 造 成 
了 天 线 的 低频 谐振 ” ， 影 响 伺服 控制 精度 ， 进 而 降低 了 天 线 跟 踪 精 度 。 
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2 伺服 控制 设计 要 求 及 算法 分 析 


2.1 伺服 控制 设计 要 求 

当天 线 系统 硬件 方案 确定 后 ， 影 响 天 线 伺服 性 能 的 最 主要 因素 是 控制 器 K(s) 的 算法 。 在 伺服 控 
制 系统 中 ， 控 制 器 K(s) 根 据 输入 设 定 值 y'(s) 与 扰动 w(s) 、 噪 声 ws) 作用 下 的 输出 值 y(s) 解 算 控制 
信号 ， 再 经 过 电机 等 驱动 装置 带动 天 线 架 转动 ， 具体 框图 如 图 1。 


扰动 w(s) 


e(s) 输出 y(s) 
C) 控制 器 K(s) 驱动 器 天 线 架 G(s) 
Wl AP) 
CO 


图 1 实际 控制 系统 框图 


Y Fig. 1 Structure of a practical control system 


输入 设 定 y'(s) 


因而 伺服 控制 设计 的 目的 是 设计 一 个 控制 器 K(s) ， 既 保证 控制 系统 的 稳定 性 、 准 确 性 以 及 快速 
Cj 性 ， 还 由 于 观测 波束 窗 ， 要 求 输出 的 超 调 量 尽 可 能 小 ， 对 一 定 的 拢 动 有 相应 的 鲁 棒 性 ， 最 终 满足 跟踪 
一 误差 的 要 求 。 由 图 1 可 得 到 系统 的 跟踪 误差 为 

C 1 G(s) K(s) G(s) 


= e(s) 2y'(s) =y(s)= (s) * (s), (1) 


1 + K(s) G(s) Y 1 tEG)G) | DeÉKGYyGG) ” 


N 由 (1) 式 可 知 ， 由 控制 器 K(s) 和 被 控 对 象 CCs) 组 成 的 开 环 传递 函数 K(s) 、C(s) 的 增益 和 带宽 需要 优 
O 化 设计 在 适当 范围 ， 使 得 系统 的 跟踪 误差 e(*) 趋 于 0， 实现 高 精度 控制 。 
CGI — 022 控制 算法 分 析 
~y 目前 应 用 或 者 研究 的 天 线 伺 服 控制 算法 (包括 对 扰动 的 处 理 ) 分 为 两 大 类 : (1) 基 于 系统 精确 模型 
的 方法 ， 比 如 线性 高 斯 二 次 型 (Linear Quadratic Gaussian, LQG) 、 内 模 控制 ( Internal. Model Control, 
IMC) 等 。 基 于 模型 建立 的 控制 律 ， 在 控制 工程 中 遇 到 很 大 的 挑战 ， 鲁 棒 性 是 首当其冲 的 问题 。 无 论 
是 基于 机 理 分 析 法 还 是 模型 辨识 法 得 到 的 系统 数学 模型 都 忽略 了 高 阶 未 建 模 动态 ， 一 旦 系统 遭遇 摄 
T O 动 , 则 控制 性 能 将 发 生变 化 。 而 基于 鲁 棒 性 提出 的 控制 算法 H_infc3 ， 其 不 确定 性 函数 在 工程 实践 中 
= 不 易 获 得 。 文 [6] 提 出 基于 LQG + GPI 的 控制 方案 可 抑制 因 大 天 线 刚度 不 足 造 成 的 拌 动 ， 但 采用 这 种 
方式 增 大 了 控制 的 调节 时 间 ， 并 不 能 给 出 最 优 的 控制 量 。(2) 非 基于 系统 模型 或 者 准确 模型 的 方法 ， 
如 传统 整数 阶 PI(D) (以 下 简写 PI(D) ) 、 分 数 阶 PID(FOPID ) 等 调节 器 。 目 前 工程 应 用 最 成 熟 的 是 
PI(D), ， 世 界 各 大 射电 望远镜 的 伺服 控制 算法 也 大 都 采用 PI 控制 。 但 是 由 于 PI 控制 器 对 采集 信号 的 
粗糙 处 理 ， 使 得 其 在 控制 精度 要 求 极 高 的 场合 达 不 到 要 求 。 虽 然 分 数 阶 PID 控制 器 阶 数 的 取 值 特点 
使 得 在 鲁 棒 性 设计 上 具有 比 传统 PID 控制 顺 更 为 灵活 的 优点 ， 系 统 不 但 可 以 满足 稳定 性 要 求 ， 同 时 
还 满足 系统 参数 鲁 棒 性 设计 的 需求 。 但 由 于 工程 实现 较 复 杂 ， 目 前 仅 限 于 理论 分 析 。 

为 实现 控制 精度 的 最 优 ， 通 过 分 析 影 响 输出 的 原因 ， 提 出 了 基于 扰动 -不 确定 性 观测 /估计 - 消 
除 '" 的 思想 ; 例如 韩 京 清 在 详细 分 析 了 控制 思想 和 传统 PID 控制 器 的 优 缺点 之 后 ， 结 合 现代 控制 中 
状态 观测 反馈 的 思想 ， 提 出 了 自 抗 扰 控 制 ( Active Disturbance Rejection Control, ADRC)。 在 自 抗 扰 控 
制 中 ,将 系统 所 受 的 外 部 扰动 和 内 部 扰动 统一 对 待 ， 借 助 扩张 状态 观测 器 (ESO) 观 测 并 补偿 扰动 ， 给 
出 控制 信和 号， 补偿 后 被 控 对 象 的 传递 函数 变 成 一 个 二 阶 积分 右 。 此 后 ， 自 抗 扰 控制 获得 了 广泛 的 研究 
和 应 用 ， 如 文 [8] 提 出 的 线性 自 抗 扰 控制 ， 直 接 将 系统 输出 值 与 观测 值 的 差 值 用 于 观测 器 的 输入 ， 形 
成 线性 扩张 状态 观测 需 (LESO ) ， 参 数 整 定 也 形成 了 一 定 的 规模 。 然 而 ， 线 性 自 抗 扰 控 制 随 着 线性 扩 
张 状态 观测 器 阶 次 的 升 高 ， 参 数 也 随 之 增加 ， 造 成 了 调试 的 困难 。 
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文 [9] 给 出 了 将 带 有 三 阶 线性 扩张 状态 观测 器 的 线性 自 抗 扰 控 制 (以 下 简称 三 阶 线性 自 抗 扰 控 制 ) 
应 用 于 高 阶 系统 的 可 能 性 以 及 参数 整定 方法 。 鉴 于 伺服 电机 的 传递 函数 不 超过 二 阶 ， 以 及 天 线 的 谐振 
频率 仅 限 制 了 控制 带宽 ， 故 考虑 使 用 三 阶 线性 自 抗 扰 控 制 。 

机 械 谐 振 的 处 理 方法 主要 有 加 速度 反馈 、 可 调 惯量 比 、 凹 口 滤波 器 (或 叫 陷 波 器 ) 和 观测 器 法 。 
加 速度 反馈 和 可 调 惯 量 比 的 实质 都 是 跟踪 、 消 除 脉动 的 驱动 力矩 。 四 口 滤 波 器 一 般 用 于 高 频 谐振 ， 借 
助 滤波 消除 谐振 ， 但 抑制 谐振 能 力 有 限 ， 且 增加 了 设计 的 复杂 性 和 难度 。 而 观测 器 法 采用 开 环 方式 估 
计 的 脉动 力矩 虽 不 完全 准确 ， 但 是 控制 效果 理想 ， 结 构 简 单 易 行 ， 故 本 文 考 虑 使 用 观测 器 法 消除 谐振 。 

通过 对 射电 望远镜 工 况 和 伺服 控制 算法 的 分 析 ， 认 为 大 口径 射电 望远镜 伺服 控制 算法 采用 三 阶 线 
性 目 抗 扰 控 制 有 较 大 的 优势 ， 文 [10-11] 也 验证 了 这 一 优势 。 
2.3 三 阶 线性 自 抗 扰 控制 与 参数 整定 

三 阶 线性 自 抗 扰 控 制 结构 框图 如 图 2。 其 中 ，y(s) 为 参考 输入 ; zi ，z,，z 为 输出 信号 y(s)、y 
(s) 微 分 值 和 系统 受到 所 有 扰动 的 估计 ; kp, ka 为 PD 控制 器 的 两 个 系数 ; bo 为 被 控 对 象 增益 的 估计 ; 
Bi, Bs. Bs 分 别 为 三 阶 线性 扩张 状态 观测 融 的 3 个 系数 ; uo(s)，u(s) 分 别 为 PD 控制 器 的 输出 和 被 
控 对 象 的 输入 信号 。 


w(s) 


被 控 对 象 


图 2 带 有 三 阶 线性 自 抗 扰 控 制 的 控制 系统 框图 
Fig.2 Structure of control system with three-order LADRC 


由 图 2 可 知 被 控 对 象 的 输入 为 
u(s) = (uo(s) -= z3) /bo, (2) 


"a uols) = kG'G) 7 2) 7 kaza (3) 
线性 扩张 状态 观测 器 的 形式 为 


sz, =z, t Bi(y(s) -21) 
SZ, = Z4 * By ys) E ; (4) 
sz, - B(y(s) — zi) 

在 估计 准确 时 ， 输 入 到 输出 的 传递 函数 G1 为 


k 
gu. e (5) 
y(s) s +hist+ k, 


z, 到 了 的 传递 函数 为 
Z3 P, 
3 3 l 6 
f P. + Bs *B,s +s? ( ) 
2E E iM SA EE TEA III, ， 即 令 : 
s? +B? +B,s +B = (s+ w,)’, 5 


Sjokkstk (sto). 
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HF, w, o, 分 别 为 线性 扩张 状态 观测 器 的 观测 带宽 和 PD 控制 器 的 控制 带宽 。 由 (7) 式 可 以 得 到 局， 
-3o,, B,73o,, B50, , k 50, k,-2v,; 改写 传递 函数 G1 为 


k, w? 
Gr 75 = 一 (8) 
s +hs+t+h, s *20, t0, 
此 时 线性 自 抗 扰 控 制 的 参数 整定 变 为 对 w.，w。 和 系统 增益 估计 b。 的 确定 。 在 时 域 分 析 中 ， 动 态 系统 


的 指标 是 调节 时 间 c, 和 超 调 量 o%， 在 阶 路 输入 下 ，(8) 式 的 超 调 量 r% 不 受 e, 影响 ， 而 根据 调节 时 
间 的 定义 ， 取 误差 为 2% 时 ，(8) 式 调节 时 间 t. 为 
t, ~ 5. 85/0, . (9) 
一 般 为 保留 一 定 的 裕 量 ，w, 应 该 大 于 上 式 求 出 值 。 
线性 扩张 状态 观测 器 的 3 个 参数 8 Ba, Bs 决定 了 观测 器 的 观测 精度 ， 经 过 极点 配置 后 3 TS 
数 是 关于 e, 的 函数 。 当 三 阶 线性 扩张 状态 观测 器 应 用 于 高 阶 系统 时 ，w。 不 能 取 太 大 ， 和 否则 导致 系统 
不 稳定 。 为 使 扩张 状态 观测 器 在 w, 取 值 较 小 的 情况 下 仍 可 以 获得 良好 的 观测 效果 ， 对 扩张 状态 观测 
器 的 参数 进行 如 下 改进 ， 即 令 : 


B, 7 3o,, B, =30,°, B. - dp,， (10) 
> rp, d 根据 被 控 对 象 的 特性 确定 。 在 调节 时 间 大 于 1 s 时 ， 一 般 可 令 : 
T 3e, > d 三 4 ， (11) 
E d 值 确定 后 ，w, 对 扩张 状态 观测 器 的 观测 速度 影响 已 经 很 小 。 一 般 可 令 ; 
T 5w, 2 w, Z 4o, . (12) 
i: 至 此 ， 线 性 自 抗 扰 控 制 除了 bo 外 的 所 有 参数 都 确定 ， 而 b, 由 调试 获得 。 


T 3 仿真 及 结果 分 析 


已 建 天 线 从 文献 中 仅 能 获得 35 m 天 线 的 传递 函数 ， 故 本 文采 用 某 35 m. 口径 天 线 对 所 选 控制 算法 
进行 仿真 验证 ， 其 方位 轴 位 置 环 传递 函数 为 
- 0, 19. 7125 + 111. 736 


G — p= A 
An U, $6.54 + 1519 + 5808? + 13705 


c eb, 0, U, 分 别 为 天 线 的 方位 输出 角 和 输入 电压 ， 利 用 MATLAB 软件 中 的 Simulink GAFA, 
建 基于 PID 和 线性 自 抗 扰 控制 器 的 控制 系统 ， 仿 真 中 取 B,/B,=6， 对 PID 与 线性 自 抗 扰 控制 中 的 PD 
和 扩张 状态 观测 器 封装 后 ， 系 统 的 仿真 框图 如 图 3。 


19.712s +111.736 
S5+6.54s4+151s3+580s2+1 370s 
方位 传递 函数 


Scope 


19.712s +111.736 
s5+6.54s4+151s3+580s?+1 370s 
方位 传递 函数 1 


图 3 基于 PID 和 线性 自 抗 扰 控制 的 仿真 框图 
Fig.3 Structure of PID & LADRC based simulation by Simulink 
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图 中 加 法 器 接 对 应 的 信和 号， 扩张 状态 观测 器 的 ina 接 测 量 噪声 。 当 输入 相同 阶 跃 信号 时 ,仿真 结 
果 如 图 4。 


-一 -- PID 


阶 跃 ”一 一 LADRC 


幅 值 


时 间 /s 
图 4 无 扰动 状态 下 阶 跃 仿真 结果 


Fig.4 Simulation result of step response without disturbance 


仿真 结果 表明 ， 不 加 干扰 时 ,线性 自 抗 扰 控制 的 调节 时 间 (2% 稳 态 误差 ) 比 PID 早 近 2s， 超 调 量 
仅 为 7% ， 为 达到 这 么 短 的 调节 时 间 ，PID 的 超 调 量 高 达 50%。 而 且 线 性 自 抗 扰 控 制 过 渡 过 程 平滑 ， 
很 好 地 解决 了 快速 性 和 超 调 之 间 的 矛盾 。 

当 加 入 扰 劲 及 噪声 时 ， 扰 劲 采 用 Davenport 功率 函数 ， 输 入 为 带 限 噪声 和 正弦 波 的 全 加 ， 用 来 模 
拟 系统 噪声 和 受到 的 波动 力矩 ， 仿 真 结果 如 图 5。 


m i 
, i 4 iM Gi W 


-1 


时 间 /s 


图 5 扰动 状态 下 阶 路 仿真 结果 


Fig.5 Simulation result of step response with disturbance 


从 图 5 可 以 看 出 ， 线 性 自 抗 扰 控 制 与 PID 都 受到 了 扰动 的 影响 。 但 线性 自 抗 扰 控 制 受 到 扰动 的 影 
响 仅 为 PID 的 25% ， 变 化 幅度 在 稳 态 误差 范围 内 ， 变 化 更 平缓 ， 表 明 控 制 性 能 优 于 PID, 

当 输 入 信号 为 带 有 一 定 初始 位 置 差 的 斜坡 信号 时 ， 用 以 模拟 天 线 跟踪 射电 源 ， 仿 真 结果 如 图 6。 
由 图 6 可 以 明显 看 到 : 线性 自 抗 扰 控制 比 PID 更 快 地 跟踪 输入 信号 ， 在 给 定 扰动 作用 下 ,线性 自 抗 扰 
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控制 的 幅 值 变化 约 为 PID 幅 值 变化 的 25%， 幅 值 变化 在 终 值 2% 的 稳 态 范围 内 ， 表 明 线 性 目 抗 扰 控 制 
不 仅 具有 良好 的 抗 扰 性 能 ， 而 且 具 有 较 高 的 跟踪 性 能 。 


5r —--- PID 参考 输入 


E 
J 
10 
ba 时 间 /s 
= 图 6 扰动 状态 下 的 斜坡 仿真 结果 
Fig.6 Simulation result of ramp response with disturbance 
4 结 d 


通过 仿真 可 以 看 到 ， 与 采用 PID 等 传统 控制 算法 相 比较 ， 线 性 自 抗 扰 控 制 带 算法 简单 易 行 ， 动 
态 特性 具有 明显 的 优越 性 。 通 过 线性 扩张 状态 观测 需 解 决 了 (1) 式 中 噪声 对 增益 和 带宽 的 限制 ， 因 而 
控制 精度 高 ， 调 节 时 间 短 ， 抗 扰 能 力 强 。 值 得 一 提 的 是 ， 本 文采 用 的 参数 整定 方法 给 出 的 是 参数 分 布 
范围 ， 具 有 更 强 的 适用 性 。 然 而 ,仿真 过 程 中 对 阵风 扰动 和 齿 隙 摩擦 的 作用 采用 了 外 加 脉动 干扰 进行 
模拟 ， 与 实际 应 用 环境 不 尽 吻合 。 虽 然 线性 自 抗 扰 控制 使 用 线性 静态 参数 扩张 状态 观测 器 显著 降低 了 
扰动 的 影响 ,但 还 不 能 做 到 无 超 调 和 无 静 差 ,该 问题 有 待 在 后 续 工 作 中 进一步 研究 。 
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Research of Control Algorithm for Antenna High Pointing Accuracy 
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Abstract: With the development of the radio astronomy and the demanding for the angular resolution, the 
pointing error of the telescope receives a very strict limitation close to 1”. In fact, it is the pulse torque that 
usually results in low frequency resonance exert to the antenna servo system, and the traditional PI ( D) 
controller could not satisfy such high pointing accuracy requirements. In this paper, considering working 
condition of radio telescope and purpose of servo control, as well as the design principle of controller, a control 
algorithm based on the Linear Active Disturbance Rejection Controller ( LADRC) is put to use. We applied this 
method to the position loop simulation of an antenna under the Simulink environment. The results prove that the 
LADRC control algorithm shows better dynamic property, higher control accuracy and stronger robustness. 


Key words: Pointing accuracy; Low frequency resonance; Servo system; LADRC ; Robustness 


